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ABSTRACT 

We investigate the evolution of the primordial gas surrounding the first massive black 
holes formed by the collapse of Population III stars at redshifts z <; 20. Carrying 
out three-dimensional hydrodynamical simulations using GADGET, we study the dy- 
namical, thermal and chemical evolution of the first relic H II regions. We also carry 
out simulations of the mergers of relic H II regions with neighboring neutral mini- 
haloes, which contain high density primordial gas that could accrete onto a Pop III 
remnant black hole. We find that there may have been a significant time delay, of 
order ~ 10 8 yr, between black hole formation and the onset of efficient accretion. The 
build-up of supermassive black holes, believed to power the z <; 6 quasars observed in 
the Sloan Digital Sky Survey, therefore faces a crucial early bottleneck. More massive 
seed black holes may thus be required, such as those formed by the direct collapse of 
a primordial gas cloud facilitated by atomic line cooling. The high optical depth to 
Lyman- Werner (LW) photons that results from the high fraction of H2 molecules that 
form in relic H II regions, combined with the continued formation of H2 inside the dy- 
namically expanding relic H II region, leads to shielding of the molecules inside these 
regions at least until a critical background LW flux of ~ 10~ 24 ergs s" 1 cm~ 2 Hz -1 
sr _1 is established. Furthermore, we find that a high fraction of HD molecules, Ahd 
<; 1CP 7 , is formed, potentially enabling the formation of Pop II. 5 stars, with masses 
of the order of ~ 10 Mq, during later stages of structure formation when the relic 
H II region gas is assembled into a sufficiently deep potential well to gravitationally 
confine the gas again. 
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1 INTRODUCTION 

What were the feedback effects from the first generation 
of stars in the Universe? The first stars, which formed in 
dark matter (DM) minihaloes of mass ~ 1O 6 M0 at redshifts 
of 2 ~ 20, were likely very massive, having characteristic 
masses of the order of ~ IOOM0 (Bromm, Coppi & Lar- 
son 1999, 2002; Abel, Bryan & Norman 2002; Nakamura & 
Umemura 2001). These massive Population III stars would 
have radiated at temperatures of ~ 10 J K (e.g. Bond, Ar- 
nett & Carr 1984), generating enough ionizing photons to 
completely ionize the minihaloes in which they were formed 
and to contribute to the reionization of the Universe (e.g. 
Barkana & Loeb 2001; Alvarez, Bromm & Shapiro 2006). 

In addition to the radiation emitted by the first stars 
during their lifetimes, those Population III stars with masses 
40 M <; M* <J 140 M or M, J> 260 M are predicted to 
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have collapsed to form black holes directly, possibly pro- 
viding the seeds for the first quasars (Madau & Rees 2001; 
Heger et al. 2003; Madau et al. 2004; Ricotti & Ostriker 
2004; Kuhlen & Madau 2005), although more massive seed 
black holes may have been formed after the epoch of the 
first stars in DM haloes with virial temperatures of Jj 10 4 
K (Bromm & Loeb 2003; Begelman, Volonteri & Rees 2006; 
Spaans & Silk 2006). 

The growth of the first black holes must have been rapid 
enough to account for the powerful quasars observed at red- 
shifts of z £ 6 (e.g. Fan et al. 2004, 2006), believed to be 
fueled by accretion onto supermassive black holes (SMBHs) 
with masses ~ 10 9 Mq (e.g. Haiman & Loeb 2001; Volon- 
teri & Rees 2005; Volonteri & Rees 2006). How such vig- 
orous accretion of matter could have taken place poses an 
important question, as it has been shown that the radia- 
tion from the first stars heats and evacuates the gas residing 
within the ~ 10 6 M minihaloes in which they are born (Ki- 
tayama et al. 2004; Whalen, Abel & Norman 2004; Alvarez 
et al. 2006). Notwithstanding some possible contribution to 
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the accreted mass from self-interacting dark matter (SIDM) 
particles (see Spergel & Steinhardt 2000; Hu et al. 2006), 
the baryonic mass around these primordial massive black 
holes (MBH) must have been efficiently replenished soon af- 
ter the birth of the black hole. In the course of hierarchical 
structure formation, this continued accretion of matter is 
naturally accomplished through mergers of the black hole's 
parent halo with its neighboring haloes (see e.g. Ricotti & 
Ostriker 2004; Kuhlen & Madau 2005; Malbon et al. 2006; 
Li et al. 2006). As recent theoretical work on the growth of 
supermassive black holes has been carried out under the as- 
sumption that Pop III seed black holes can begin accreting 
at the Eddingtion limit very soon after their formation (e.g. 
Li et al. 2006; Malbon et al. 2006), it stands as an important 
task to determine in which environments this might actually 
be possible. 

The effects of the radiation from the first stars may have 
enhanced subsequent star formation through the production 
of H2 inside relic H II regions, as well as in partially ionized 
shells just ahead of ionization fronts, as proposed by Ricotti, 
Gnedin & Shull (2001, 2002; sec also Ahn & Shapiro 2006). 
The former possibility has recently received considerable at- 
tention (Oh & Haiman 2003; O'Shea et al. 2005; Nagakura 
& Omukai 2005; Johnson & Bromm 2006). O'Shea et al. 
(2005) have reported that second-generation star formation 
could have occurred in the ionized minihaloes neighboring 
the first stars, owing to the formation of H2 molecules in the 
recombining primordial gas (see Shapiro & Kang 1987; Fer- 
rara 1998; Ricotti, Gnedin & Shull 2001, 2002). However, 
Alvarez et al. (2006) find that neighboring minihaloes are 
self-shielded to the ionizing radiation of the first stars, and 
thus that star formation in neighboring minihaloes may not, 
in fact, have been significantly enhanced. Also, it has been 
shown that the activation of cooling by deuterium hydride 
(HD) molecules inside relic H II regions may provide an av- 
enue for the formation of Population II. 5 (Pop II. 5) stars, 
with masses of the order of 10 Mq and formed from strongly 
ionized primordial gas (Mackey, Bromm & Hernquist 2003; 
Johnson & Bromm 2006; see also Nagakura & Omukai 2005). 
However, it remains to fully elucidate the formation process 
of Pop II. 5 stars within the first relic H II regions if, indeed, 
the radiation from the first stars evacuates the gas contained 
in their parent haloes and, yet, does not substantially ionize 
the gas in its neighboring minihaloes. 

Here we present the results of three-dimensional numer- 
ical simulations of the recombination of the first relic H II 
regions, investigating the possibility of Pop II. 5 star forma- 
tion in such regions. We assume, for this case, the first star 
to have a mass of 100 M and to collapse directly to a black 
hole. Additionally, we simulate the merger of this parent 
halo with a neighboring neutral minihalo which has not yet 
experienced star formation, in order to determine the neces- 
sary conditions for the black hole to begin accreting gas at 
the Eddington limit, and so to grow to a mass of 10 9 M by 
a redshift of z ~ 6. In future work, we will study the feed- 
back from a pair-instability supernova, the other possible 
fate predicted for single primordial stars with masses 100 
Mq (e.g. Rakavy & Shaviv 1967; Bond, Arnett & Carr 1984; 
Heger et al. 2003). The details of our numerical methodol- 
ogy are given in Section 2. The results of our simulations 
of the recombination of the relic H II region appear in Sec- 
tion 3. Our results from the simulation of the merging of 
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Figure 1. Comparison of results from our one-zone model 
(see Johnson & Bromm 2006) with the results from a three- 
dimensional numerical simulation with GADGET, for the case 
of a collapsing spherical minihalo of uniform density. The solid 
lines show the results of the one-zone model, while the triangu- 
lar symbols show output from the three-dimensional simulation. 
There is clearly good agreement between the two calculations for 
the evolution of the gas density (top- left panel), as well as for the 
H2, HD, and the free electron abundances (top-right, bottom-left, 
and bottom-right panels, respectively). 

the relic H II region with a pre-collapse halo are presented 
in Section 4, while the implications of these results for the 
growth of the remnant black hole appear in Section 5. Fi- 
nally, in Section 6 we summarize our results and present our 
conclusions. 



2 METHODOLOGY 
2.1 Chemical network 

We employ the parallel version of the GADGET code for 
our three-dimensional numerical simulations. This code in- 
cludes a tree, hierarchical gravity solver combined with the 
smoothed particle hydrodynamics (SPH) method for track- 
ing the evolution of the gas (Springel, Yoshida & White 
2001). Along with H 2 , H 2 + , H, H", H+, e", He, He+, and 
He ++ , we have included the five deuterium species D, D + , 
D~, HD and HD~, using the same chemical network as in 
Johnson & Bromm (2006). 

As a test of the reliability of the chemical network in- 
corporated into GADGET, we have simulated the idealized 
case of the homologous collapse of a ~ 10 M spherical 
cloud with an initial uniform density of m ~ 2 cm' 3 . The 
chemical species are initialized with their primordial abun- 
dances, as given in Galli & Palla (1998). The temperature of 
the gas, initially a uniform 200 K, was chosen so that there 
would be little pressure support of the gas against gravita- 
tional collapse, allowing the cloud to collapse essentially in 
free-fall. We then followed the thermal and chemical evo- 
lution of the gas near the center of the sphere, where the 
density profile is approximately flat. We compare our simu- 
lation results with those obtained from the one-zone model, 
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Figure 2. The properties of the primordial gas in the minihalo 
identified to host the first star, as functions of distance from the 
center. The values shown for the central density, temperature, 
c~ fraction, and H2 fraction are very close to the canonical val- 
ues generally found in simulations of Pop III star formation (e.g. 
Bromm & Larson 2004). 



employed in previous work (see Johnson & Bromm 2006). In 
this one-zone calculation a cloud of uniform density collapses 
homologously under its own gravity, its density evolving ac- 
cording to 



dn 
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(1) 



where n is the density at time t after the onset of the col- 
lapse, 771h is the mass of the hydrogen atom, and n is the 
mean molecular weight. Fig. 1 shows a comparison of the 
evolution of the density and of the abundances of H2, HD, 
and free electrons found from our GADGET simulation with 
that found from the calculation using our one-zone model. 
The agreement is very good, giving us confidence in the ac- 
curacy of our chemical network. 



2.2 First ionizing source 

The initial conditions for our three-dimensional SPH cal- 
culation are given by a cosmological simulation of high-z 
structure formation that evolves both the dark matter and 
baryonic components, initialized according to the ACDM 
model at z — 100. In carrying out the cosmological simula- 
tion used in this study, we adopt the same parameters as in 
earlier work (Bromm, Yoshida & Hernquist 2003). We thus 
use a periodic box of size L = 100 h~ x kpc comoving and a 
number of particles Nom = Nsph = 128 3 . The SPH particle 
mass here is ~ 8 Mq. 

We identified the first SPH particle to achieve a density 
above 10 4 ' 5 cm~ 3 within our cosmological box at z ~ 19.5, 
finding it at the center of a minihalo with a total mass of 
~ 10 6 Mq. Fig. 2 presents the properties of the primordial 
gas as a function of distance from the minihalo center. The 
gas temperature rises as particles are adiabatically heated as 
they fall into the potential well of the halo, and then drops 



Figure 3. The density, temperature and radial velocity of pri- 
mordial gas ionized and heated by radiation from the first star, 
as a function of distance from the star, for two representative 
stellar masses. Here we show the situation after 3 Myr for the 
case of the 100 Mq star and after 2 Myr for the case of the 200 
Mq star. The 100 Mq star will likely collapse directly to form a 
black hole after this time, while the 200 Mq is likely to explode 
as a pair-instability supernova. In the present work, we track the 
evolution of the ionized gas as it recombincs and cools in the case 
of the 100 Mq black hole-forming star. 



nearer the center of the halo where the H2 fraction rises to 
~ 10~ 3,4 and molecular cooling can thus efficiently cool the 
gas to ~ 200 K (e.g. Bromm & Larson 2004). 

Having identified the location of the first star, we placed 
a point source of ionizing radiation at that location in our 
cosmological box. This was effected by including the fol- 
lowing heating rates and ionization rate coefficients in our 
calculations of the thermal and chemical evolution of the 
gas: 
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where r is the distance from the star in pc, and the subscripts 
denote the chemical species subject to photoionization and 
photoheating. These heating rates and ionization coefficients 
are derived from the models given in Schaerer (2002) for the 
case of a ~ 100 Mq Pop III star, assuming the stars emit 
a blackbody spectrum (see also Bromm, Kudritzki & Loeb 
2001). We also carried out a simulation of the ionization 
and heating of the gas assuming a stellar mass of 200 Mq , 
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and the resulting density, temperature, and radial velocity 
profiles around the central ionizing source, for each case, are 
shown in Fig. 3. The profiles show the situation after 3 Myr 
and 2 Myr of photoheating and photoionization, for the f 00 
Mq and 200 Mq cases, respectively (see e.g. Schaerer 2002). 
The higher effective temperature and luminosity of the 200 
Mq star results in both a harder spectrum and more ionizing 
photons, and so in a higher heating rate of the surrounding 
gas from photoionization. Thus, as can be seen in Fig. 3, 
the temperature of the gas at a given distance from the 
central source is at least several ~ 10 3 K higher for the 
case of the 200 M star than for that of the 100 M star. 
Also, due to the shorter lifetime of the 200 Mq star, the 
shock that arises from the steep temperature and density 
gradients encountered during the photoheating of the gas 
has not moved as far out from the central star for the 200 
M case, although the shock velocity is higher in this case. 

Although we neglect the detailed radiative transfer of 
the ionizing photons here, we succeed at reproducing the 
basic features of the temperature and density profiles of 
the gas around the ionizing source that have been found 
in previous radiative transfer calculations (Kitayama et al. 
2004; Whalen et al. 2004; Alvarez et al. 2006). Also, we heat 
and ionize the gas only within 500 pc of the central source, 
which is roughly consistent with an inhomogeneously ion- 
ized region around the first star of order ~ 1 kpc, as found 
by Alvarez et al. (2006), without impinging on neighbor- 
ing minihaloes in our cosmological box. We require that our 
H II region not encompass any neighboring minihaloes be- 
cause these can be self-shielded to ionizing radiation even if 
they reside inside the H II region, and so we cannot accu- 
rately follow the chemical evolution of the gas inside those 
minihaloes while the central ionizing source is on. Since we 
do not explicitly follow the propagation of the ionization 
front with time, we also do not resolve the time-dependent 
effects on the chemistry and thermal evolution of the gas 
that can give rise to, for instance, the formation of shells of 
H2 molecules just outside the I-front (see Ricotti, Gnedin 
& Shull 2001, 2002). Despite the inability of our method to 
capture the detailed structure of the H II region, we do ex- 
pect that we can capture the essential chemical and thermal 
evolution of the relic H II region as a whole in our simula- 
tions, as it expands, cools, and recombines after we remove 
the central source from the calculation. 

We used the version of GADGET which integrates the 
entropy equation for our photoionization calculation, as we 
also do later for our simulations of the relic H II region. As 
opposed to the integration of the energy equation, this for- 
mulation of GADGET conserves both energy and entropy 
and is much more successful at resolving the thermal evo- 
lution of gas that experiences shocks or strong local energy 
injection (Springel & Hernquist 2002). 

2.3 Recombination and molecule formation 

With the formation of a black hole by direct collapse of the 
100 Mq Pop III star, the relic H II region left behind begins 
to recombine and cool. We implement this by simply setting 
the photoionization coefficients and heating rates to zero, for 
the case of the 100 Mq central star. Thus, the temperature, 
density, and radial velocity profiles shown in Fig. 3 are the 
initial conditions for our simulation of the relic H II region. 



We follow the thermal and chemical evolution of the relic 
H II region for ~ 100 Myr, considering in particular the 
production of molecules and the cooling of the primordial 
gas. 



2. 3. 1 Photodissociation of molecules 

Molecular hydrogen can easily be dissociated by absorption 
of Lyman- Werner (LW) photons with energies between 11.2 
and 13.6 eV (e.g. Haiman, Rees & Loeb 1997; Bromm & 
Larson 2004, and references therein). Although it is well- 
established that an external LW background could be pro- 
duced by stars born in neighboring minihaloes (e.g. Haiman, 
Rees & Loeb 1997; Haiman, Abel & Rees 2000; Ciardi et al. 
2000; Machacek, Bryan & Abel 2001), here we assume that 
at the death of our first star a negligible UV background 
has been established by emission from stars elsewhere in 
the Universe. Thus, to evaluate the effect that photodisso- 
ciation has on the molecule fraction in the first relic H II 
region, we consider as the only source of dissociating radia- 
tion two-photon emission from the 2 1 S — > 1 1 S transition in 
recombining helium atoms from within the relic H II region 
itself (Johnson & Bromm 2006). Given the much larger Ein- 
stein A coefficient for two-photon emission from 2 1 S than 
from 2 3 S, 51.3 s" 1 for the 2 1 S -> 1 1 S transition versus 2.2 
x 10~ 5 s _1 for the 2 3 S -> 1 1 S transition, this should be a 
sound approximation (Mathis 1957; Osterbrock & Ferland 
2006). 

To include a prescription for the photodissociation rates 
of H2 and HD in our code we assume for simplicity that the 
relic H II region is spherically symmetric and that it is opti- 
cally thin to the LW photons. With this latter assumption, 
we obtain an upper limit for the dissociation rate, as the 
molecule fraction can approach ~ 10~ 3 in relic primordial 
H II regions, which may lead to an appreciable optical depth 
to LW photons (see e.g. Ricotti, Gnedin & Shull 2001; Oh & 
Haiman 2003; Kuhlen & Madau 2005; O'Shea et al. 2005). 
We estimate the total number of He recombinations, He + + 
c~ — > He + hv, per second within the H II region that lead 
to population of the 2 X S state, Q 2 1 S) according to 



Q 



2 1 S 



E QBn c nHciimspH 
3/imHn 



(8) 



where the sum is over all SPH particles in the H II region. 
Here n c is the number density of free electrons, nn e 11 is the 
number density of He II, n is the total number density, mspH 
is the mass per SPH particle, \x is the mean molecular weight, 
and ob is the Case B total He recombination coefficient to 
singlet states. We take it that 1/3 of the recombinations 
to the singlet levels of He I result ultimately in population of 
the 2 1 S state, which is accounted for by the factor of 1/3 in 
the above formula (see Pottasch 1961; Osterbrock & Ferland 
2006). For the LW flux at the edge of the H II region, we 
find 

; Q2IS 



Jlw ~ 10 



4irR' 2 ' 



(9) 



where R is the radius (in cm) of the He II region, Jlw is the 
LW flux in units of 10~ 21 erg s _1 cm~ 2 Hz -1 sr _1 , and we 
have conservatively estimated the probability of LW photon 
emission per two-photon transition 2 1 S — > 1 X S as ^ 0.4 (see 
Osterbrock & Ferland 2006) . We compute the timescale for 
the photodissociation of H2 and HD as 
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where we have taken R ~ 500 pc and used idiss ~ 10 yr 
(^lw/10 -4 ) -1 (see Oh & Haiman 2003; Johnson & Bromm 
2006). Inverting this timescale, we find a typical rate for the 
dissociation of H2 and HD of 

10 " 16 S_1 (l^) • (") 

which is included in the calculation of the molecule fraction 
in our simulations. 

To give an estimate of the importance of the dissociation 
of molecules due to two-photon emission, we calculate t diss 
for the simplified case of a spherical recombining H II region 
of uniform density. In this case, we have the total number 
of recombinations to He I which ultimately result in the 2 1 S 
state given by 

4ir 3 

^ 10 » B -i(_j?_y (im^l) . ( i2) 

\500pc J V cm- 6 J y ' 

Using equation (10), we find 

t diss ~ 10 yr - j , (13) 

which clearly shows that at high densities and at times when 
He recombination is still ongoing, the photodissociation of 
any molecules that have formed may be important, even 
in the absence of an externally generated LW background. 
However, for our case of a relic H II region surrounding a 
minihalo, since it is predominantly at lower temperatures (^ 
5,000 K) that molecules are formed, at which times much of 
the He II has already recombined, and since our H II region 
is expanding to lower densities with time, we expect that 
dissociation of molecules due to two-photon emission will be 
unimportant, at least during the later evolutionary stages. 

Molecular hydrogen could also be dissociated by radia- 
tion generated during the accretion of gas onto the remnant 
black hole. However, as we show in Section 5, we find that 
the accretion rate onto the black hole is low, comparable to 
that found by O'Shea et al. (2005), for at least a few 10 Myr 
after the collapse of the central ionizing star. O'Shea et al. 
estimate that this accretion rate results in a photodissocia- 
tion rate of H2 which is at least an order of magnitude below 
the formation rate of H2 in the relic H II region. We thus ne- 
glect the possible effects of photodissociating radiation due 
to accretion onto the remnant MBH. 

2.4 Merging minihaloes 

In order to determine which, and how quickly, neighbor- 
ing minihaloes will collapse following the formation of a 
Pop III remnant black hole, we consider the situation in 
which the relic H II region surrounding the remnant black 
hole merges with a neighboring DM halo and its accompany- 
ing, un-ionized and dense gas component. Since we are sim- 
ulating the evolution of the relic H II region and the infalling 
neutral minihalo after the death of the Pop III star, in these 
merger simulations we do not include any photoheating or 
photodissociating terms in our calculations of the thermal 



and chemical history of the gas. We initiate the merger by 
placing the spherical relic H II region, immediately follow- 
ing the collapse of the central star, with properties shown in 
Fig. 3, adjacent to a second spherical region of radius 500 pc 
at the center of which is a minihalo which is still neutral, 
not yet having hosted star formation. The region containing 
the neutral minihalo is selected and cut out from elsewhere 
in the same cosmological box, and is then placed adjacent to 
the relic H II region in a new, otherwise empty simulation 
box. The initial separation between the centers of the two 
minihalos is 1 kpc proper for all merger simulations carried 
out here. 

The reason why we must choose these initial conditions, 
and why we may not simply continue running our cosmolog- 
ical simulation of the relic H II region and wait for a merger 
to occur, is that we are limited by the size of the cosmologi- 
cal box. Our box size is ~ 100fe _1 kpc, which is too small to 
contain the large wavelength density modes that drive the 
mergers of minihaloes. Thus, we carry out the mergers in 
an empty box, imparting a relative velocity of 8 km s _1 to 
the merging halos, comparable to the virial velocity of the 
resulting system, setting them on trajectories for a direct 
collision. 

We assume that the molecules in the infalling, pre- 
collapse halo have been destroyed by the H2 photodissociat- 
ing LW flux from the nearby Pop III star which has formed 
the relic H II region. We carry out simulations in which 
the neighboring halo has peak gas densities of ~ 0.1 cm -3 , 
1 cm~ 3 , 10 cm~ 3 , and 100 cm~ 3 , corresponding to different 
degrees of pre-collapse. In the case of the 0.1 cm -3 peak den- 
sity neighboring halo, the free-fall time will be comparable 
to the timescale for the completion of the merger and so we 
can neglect the possibility of this minihalo collapsing to form 
a star before the completion of the merger. However, given 
that we assume that there are no molecules in the neighbor- 
ing halo at the outset of the merger, we expect that collapse 
will be delayed even for the higher density cases, as cool- 
ing will be suppressed until molecules have reformed (see 
Mesinger et al. 2006). 

Although the initial conditions for these simulations are 
idealized, as we have not followed the merger of the haloes in 
a fully cosmological context but instead in a box containing 
only the relic H II region and the infalling minihalo, we 
are able to discern the crucial aspects in the thermal and 
chemical evolution of the gas in the vicinity of the remnant 
Pop III black hole. We do note, however, that additional 
effects that we do not consider here, such as the formation 
of an H2 shell and the driving of a shock through a partially 
ionized minihalo, can become important for cases in which 
significant portions of the infalling minihalo are ionized (see 
Ahn & Shapiro 2006). A more sophisticated treatment of 
radiative transfer will be required to more accurately follow 
the evolution of the gas (see Ahn & Shapiro 2006; Susa & 
Umemura 2006) . 



3 EVOLUTION OF PRIMORDIAL GAS IN 
RELIC H II REGIONS 

The results of our recombination simulation are presented 
in Fig. 4 at three representative times after the death of the 
central star. As can be seen in the panels showing the tem- 
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perature as a function of the density, the gas initially cools 
largely by adiabatic expansion, as at temperatures below 
~ 10 4 K the gas closely follows the adiabatic relation T oc 
n 2,/3 , delineated in the panels on the right-hand side. As the 
gas recombines, the cooling rate due to collisional excitation 
of the newly-formed hydrogen atoms is enhanced and the 
temperature of the gas drops to ~ 10 4 K, at which point the 
cooling rate decreases when molecular hydrogen becomes the 
main coolant, aside from adiabatic cooling, which continues 
as the gas expands into the intergalactic medium (IGM). 

The HD fraction in the highest density regions of the re- 
combining gas increases to Xhd ~ 10~ 7 after ~ 10 Myr since 
the death of the central star, and to Xhd ~ 10~ 6 after 100 
Myr. Thus, the HD fraction quickly rises above the critical 
value of XHD.crit ~ for efficient cooling of primordial 

gas in local thermodynamic equilibrium (LTE) to the tem- 
perature of the cosmic microwave background (CMB) (see 
Johnson & Bromm 2006) . Because LTE can only be achieved 
at much higher densities than those that persist in the relic 
H II region we consider here, radiative cooling to the CMB 
floor would only be a viable possibility if densities are some- 
how increased to the point that LTE can be established. 
This would happen if the gas is at some point incorporated 
into a larger DM halo and becomes gravitationally bound 
once more. That the HD fraction is so high, however, means 
that the potential for Pop II. 5 star formation does exist, in 
principle, if the gas becomes gravitationally bound and thus 
available for star formation (see also Nagakura & Omukai 
2005). 

We evaluate the recombination time after 100 Myr by 
taking the density of H + and of free electrons to be n H + ~ 
n e ~ 10~ 4 cm~ 3 , as can be seen from the lower-left panel 
of Fig. 4. Then, assuming at these low densities a Case A 
recombination coefficient of a a ~ 6 x 10 -13 cm 3 s _1 , we 
obtain a recombination time of t ICC ~ 500 Myr. This is more 
than twice the Hubble time at these redshifts, and suggests 
that the free electron fraction left over from the ionization 
caused by the first stars may have remained an important 
catalyst for molecule formation even after hundreds of mil- 
lions of years since the death of the central star (see e.g. 
Shapiro & Kang 1987; Yamada & Nishi 1998; O'Shea et al. 
2005; Nagakura & Omukai 2005). 

The optical depth to LW photons becomes unity for H2 
column densities of ~ 10 14 cm~ 2 (e.g. Draine & Bertoldi 
1996; Osterbrock & Ferland 2006), and for our relic H II 
region we estimate 



TLW 



nnXn 2 R 
10 14 cm- 2 



(14) 



where Xn 2 is the molecule fraction, R is the radius of the 
relic H II region, and nu is the number density of hydrogen 
nuclei. We find that the molecule fraction approaches Xu 2 ~ 
10~ 3 and that the number density becomes riu ~ 10 _3 cm~ 3 . 
Taking the radius of the H II region to be ~ 1 kpc, we find 
that the optical depth to LW photons becomes of the order 
of tlw ~ 10. If the density of star-forming minihaloes is thus 
high enough, an appreciable suppression of the background 
LW flux may result from the high molecule fraction which 
arises in relic H II regions. This may provide an important 
degree of shielding from molecule-dissociating radiation and 
may lead to a higher overall efficiency of star-formation in 
minihaloes with virial temperatures <, 10 4 K (see Ricotti, 



Gnedin & Shull 2001; Machacek, Bryan & Abel 2001, 2003; 
Oh & Haiman 2002). 

We carried out simulations both with and without the 
photodissociating two-photon emission from recombining He 
II included, and we found that this photodissociating radi- 
ation had little effect on the molecule abundances. To es- 
timate the level of LW background radiation necessary to 
efficiently photodissociate H2 molecules, we evaluate the H2 
formation timcscale at ~ 10 Myr after the death of the cen- 
tral star. Taking representative values for the temperature, 
number density, and abundances of the chemical species af- 
ter 10 Myr of recombination, we find a formation timescale 
for H2 of tform ~ 10 7 yr. The continued formation of H2 
is driven by the high abundances of H, H~, and e~, which 
are the reactants in the following reaction sequence that is 
the most important for the production of H2 (e.g. Kang & 
Shapiro 1992): 



H + e~ 



H + H" 



H - + hv 



H 2 +e" 



(15) 



(16) 



Equating the formation timescale for H2 with the dissocia- 
tion timescale for H2, given by i d i ss ~ 10 8 yr (Jlw/10 -4 ) -1 , 
we find a critical value for the background LW flux, below 
which molecules in the relic H II region will not be pho- 
todissociated efficiently, of the order of JLW.crit ~ 10 (see 
also Oh & Haiman 2003). This is comparable to the back- 
ground LW flux that is expected to have been established by 
the first generations of stars at redshifts z i; 15 (e.g. Greif 
& Bromm 2006). Taking into account the optical depth to 
LW photons of order tlw ~ 10, we find that the most heav- 
ily self-shielded molecules in the center of relic H II region 
could only be dissociated by a background LW flux, emanat- 
ing from outside the relic H II region, at least of the order 
of Jlw ~ 10. This value would, however, decrease with time 
if the molecules nearer to the periphery of the relic H II re- 
gion were dissociated by the external LW background flux, 
and so become unavailable for shielding the inner molecules 
from the dissociating radiation. We note also that if the gas 
in the relic H II region has large velocity gradients, owing 
to turbulence that we do not resolve in these simulations, 
then the optical depth to LW photons may be lower than 
the value we find here by a factor of a few (see Draine & 
Bertoldi 1996; Osterbrock & Ferland 2006). This would not, 
however, effect the value that we find for JLw.crit, as this is 
independent of the optical depth to LW photons. 

The values that we find for JLw.crit and for tlw suggest 
that the enhanced fractions of H2 and HD inside relic H II 
regions could have persisted down to at least redshifts of 
z ~ 15, and so, in principle, would have been available for 
star formation at least down to these redshifts (see also Ri- 
cotti, Gnedin & Shull 2002). The high HD abundance, which 
becomes at least an order of magnitude above the critical 
abundance for cooling the primordial gas to the CMB tem- 
perature floor, could thus have led to Pop II. 5 star formation 
inside the first relic H II regions, if these regions became in- 
corporated into more massive DM haloes that could gravita- 
tionally bind, and so increase the density of, the recombining 
primordial gas. 
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Figure 4. The evolution of the relic H II region. From left to right, the panels show the free electron fraction, the H2 fraction, the HD 
fraction, and the temperature as functions of density at ~ 1 Myr (top row), ~ 10 Myr (middle row), and ~ 100 Myr (bottom row) 
after the collapse of the central star to a black hole. The long-dashed line in the rightmost panels denotes the temperature of the cosmic 
background radiation, 7cmb> while the short-dashed line denotes the line T oc n 2 / 3 , along which gas evolves adiabatically. Here, we plot 
only the SPH particles that were subjected to the photoionizing radiation of the central star, that is, particles within ~ 500 pc of the 
central star. 



4 EVOLUTION OF PRIMORDIAL GAS IN 
MERGING MINIHALOES 

To discern the conditions under which the Pop III remnant 
black hole could efficiently accrete the dense, cold gas sup- 
plied by a neighboring minihalo, we have tracked the evo- 
lution of the gas in such merging systems with a range of 
initial peak densities of the gas within the infalling neu- 
tral minihalo. Again, here we have assumed that the LW 
flux from the now-collapsed Pop III star has destroyed all 
of the molecules inside the neutral infalling halo, although 
molecules can reform during the merger, since there is no 
longer a LW flux from the original Pop III star. The initial 



peak densities of the halos we follow are 0.1, 1, 10, and 10 2 
cm~ 3 . Fig. 5 shows the evolution of the density structure in 
the merger between the relic H II region and a halo with a 
peak density of 10 cm~ 3 , as a representative case. 

Fig. 6 shows the time evolution of a merger of a relic 
H II region, in which the central star has just collapsed to 
form a black hole, with a neutral pre-collapse minihalo which 
has a peak density of 0.1 cm~ 3 at the time of the formation 
of the black hole, at four representative times. As in the case 
of the recombination of the relic H II region evolved in our 
100/t -1 kpc cosmological box, which did not experience a 
merger, the relic H II region gas expands and cools largely 
adiabatically. This expansion is evident in Fig. 5 as well. 
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Figure 5. The merging of the relic H II region with a neutral neighboring minihalo with an initial central density of ~ 10 cm -3 . The 
H II region is on the right-hand side of the top-left panel at the beginning of the merger, while the prc-collapsc halo with which it merges 
is on the left-hand side. The location of the remnant black hole, initially at the center of the relic H II region, is shown by the black 
square in each panel. Here, we assume that the black hole has a ballistic trajectory, with a constant velocity of 4 km s — 1 to the left, 
equal to the initial velocity of the relic H II region. The highest density gas is shown in white and the lower density gas is shown in blue. 
The merger is shown at 1 Myr (top-left), 10 Myr (top-right), 50 Myr (bottom- left), and 100 Myr (bottom-right) after the death of the 
central star. The halo shown on the left-hand side collapses to a density of ~ 10 3 cm -3 during the merger (see Fig. 7). 



Also, as can be seen from the temperature rise in the pre- 
collapse halo, shown in black in Fig. 6, the expansion of 
the relic H II region, combined with the 8 km s" 1 relative 
velocity of the infalling halo with respect to the relic H II 
region, shock-heats the neutral gas within the infalling halo, 
contributing to the suppression of the density in the pre- 
collapse halo. The gas in the pre-collapse halo, furthermore, 
does not reform H2 molecules efficiently, owing to the low 
density of the gas in this halo. The H2 fraction approaches 
only ~ 10" 6 after 100 Myr since the death of the central 
star. Thus cooling of the gas is inefficient, and the highest 
densities achieved in this merger are n ^ 10 -0 ' 5 cm -3 . 



We can expect the gas within a minihalo that merges 
with a relic H II region to be dispersed and have its density 
suppressed whenever the ram pressure of the expanding relic 
H II region, given by P ram ~ nu nmnv 2 , is higher than the 
pressure of the neutral gas within the infalling halo, given by 
P gas = ngasfcBT, where v is the expansion velocity of the relic 
H II region, iih n is the density in the relic H II region, n gas 
is the density of the gas within the neutral minihalo that is 
merging with the relic H II region, and T is the temperature 
of this gas. Taking a fiducial value of v ~ 10 km s _1 for the 
expansion velocity, as can be seen in Fig. 3, and taking nn 11 
~ 10 -2 cm -3 as a typical value for the density of the relic 
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Figure 6. The evolution of the merger of the relic H II region 
with a neutral neighboring minihalo. The H II region gas is in 
orange, while the gas from the pre-collapse minihalo is in black. 
The central density of the pre-collapse minihalo is initially ^0.1 
cm -3 , roughly 10 2 times the background IGM density. We follow 
the same convention for the lines delineating the CMB tempera- 
ture and adiabatic evolution as in Fig. 4. The criterion given by 
equation (17) for neutral gas retention is not satisfied in this case, 
and the low density gas in the infalling halo is shock-heated and 
remains at low densities during the merger. Furthermore, the H2 
fraction stays below 10 — 6 for the entire 100 Myr duration of the 
merger, further preventing the gas from cooling and collapsing to 
higher densities. 

H II region, as can be seen in Fig. 4, we find the following 
condition for the retention of the neutral gas in the potential 
well of the infalling halo: 

n gas r^ 10 2 Kcm- 3 . (17) 

The condition in equation (17) can be satisfied for primor- 
dial gas collapsing in a minihalo with a mass of ~ 10 6 Mq, 
provided that densities of at least n ~ 10 -0 ' 5 cm -3 have 
been reached prior to the merger [see Fig. 10 in Bromm et 
al. (2002)]. 

However, if the gas in the infalling minihalo collapses 
to form a star before the completion of the merger, then 
the gas in the minihalo will be heated and expand to lower 
densities, as shown in Fig. 3 (see also Abel, Wise & Bryan 
2006). This final collapse to form a star will occur on the 
order of the free-fall time, ts oc n _1//2 , if the free-fall time 
is longer than the cooling time (e.g. Tegmark et al. 1997; 
Ciardi & Ferrara 2005). If the gas cannot cool efficiently, 
however, this collapse will be delayed. This may indeed be 
the case in regions near the first stars where the LW flux 
generated during their ^ 3 Myr lives could destroy the H2 
molecules inside the pre-collapse haloes, depriving them of 
the coolants that allow for star formation (e.g. Yoshida et 
al. 2003). 

To find the highest densities that could be achieved dur- 
ing a merger of a neutral minihalo with a relic H II region, we 
compare the timescale for completion of the merger, £ m er g e, 
with the timescale for the collapse of the neutral minihalo, 
tcoiiapse, the latter found from our simulations of mergers 



Figure 7. Same as in Fig. 6, except that the neutral gas has an 
initial peak density of ~ 10 cm -3 . The criterion in equation (17) 
for neutral gas retention is satisfied here, and the densest gas, 
in the center of the infalling minihalo, retains its high density 
despite the outer layers of gas being shock-heated in the merger. 
Also, the H2 fraction approaches 10 — 3 after 100 Myr at these 
higher densities, so that the gas in the merging halo cools and 
collapses efficiently to a density of ~ 10 3 cm -3 . As the labels 
indicate, the four panels correspond to the same times since the 
death of the central star as the four panels in Fig. 5, which shows 
the evolution of the projected gas density. 

involving minihaloes with peak densities of 1, 10, and 10 2 
cm -3 at the time of the formation of the black hole and the 
cessation of the radiation from the original Pop III star. For 
each of these initial densities, the criterion given by equation 
(17) for neutral gas retention is satisfied. Cases for which 
tmerge ^5 tcoiiapso will give rise to mergers resulting in the 
highest densities of gas that can be accreted onto the black 
hole, as it is these minihalos that will merge completely with 
the black hole before collapsing to form a second star. We 
define the merger timescale as 

f° d r 

Emerge — / " " j~ , (18) 

7,. MBH (2GM halo (i - ^))i 

where r MBH is the initial distance between the centers of the 
merging haloes at the time of the formation of the massive 
black hole in the relic H II region, and Mu&ia ~ 1O 6 M0. 
We note, however, that this is only an approximation to 
the actual time that would be required for a merger to take 
place, as this formula assumes that the merging minihaloes 
start at rest with respect to each other at the time of the 
collapse of the first star, and this will not be the case in 
general (see e.g. Abel, Wise & Bryan 2006). 

Fig. 7 shows the time evolution for a merger with a neu- 
tral peak density of 10 cm -3 , at four representative times, 
just as in Fig. 6. The gas in the pre-collapse halo, in this 
case, does reform H2 molecules efficiently, owing to the high 
density of the gas in this halo. The H2 fraction approaches 
~ 10~ 3 after 100 Myr. Thus, gas cooling is efficient, and the 
density reaches n ~ 10 3 cm -3 after ~ 100 Myr. Fig. 5 shows 
the evolution of the gas density structure in this merger, at 
the four representative times which are also shown in Fig. 7. 
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Figure 8. Requirements for the efficient accretion of gas onto a 
Pop HI remnant black hole. The two horizontal lines show the 
time it takes for infalling minihaloes with central densities of 10 
and 100 cm -3 to collapse and form stars. The merger timescale, 
tmerge, defined in equation (18), is the time it takes for the rem- 
nant black hole to merge with the infalling minihalo, and is a 
function of rjyiBH, the distance between the black hole and the 
center of the infalling minihalo at the time of the formation of the 
black hole. For the black hole to efficiently accrete gas at near the 
Eddington rate, it must merge with the infalling halo before this 
halo collapses to form a star. Thus, efficient accretion onto the 
black hole is only possible if t m0 rgc < tcollapsc- 



For the case of a merger with an initial peak density of 
1 cm" 3 , we find that the H2 fraction becomes only of the 
order of 10 -5 , and the density in this minihalo thus does 
not increase beyond ~ 1 cm" 3 within 100 Myr, since molec- 
ular cooling is suppressed. We also find that for the case of a 
merging minihalo with an initial density of 100 cm -3 the H2 
fraction becomes ~ 10" 3 and the halo collapses after ~ 60 
Myr. As expected, for higher initial peak densities in the in- 
falling haloes, the timescales for the collapse of these haloes 
become shorter, both because H2 molecules are reformed 
more quickly and because the free-fall time is shorter for 
denser haloes (see also Mesinger et al. 2006). 

Fig. 8 shows the requirements for infalling minihaloes 
to collapse to high densities on a timescale that is longer 
than the timescale for their merging with a black hole inside 
a relic H II region. These infalling minihaloes, which have 
Emerge & tcollapsc, are those that will contain high density 
gas (<; 10 2 cm" 3 ) at the completion of the merger, which 
can lead to efficient accretion onto the black hole. As can 
be seen in Fig. 8, at the time of the death of the original 
Pop III star, these minihaloes have densities <; 10 cm" 3 and 
lie within tmbh ~ 340 pc of the halo hosting the first Pop III 
star. Higher densities than this could be achieved, within a 
merger timescale of the order of ~ 100 Myr and without 
forming a second star, in mergers with minihaloes that have 
initial densities of the order of 10 2 cm -3 only if the timescale 
for completion of the merger is t m cr g c 60 Myr, the time 
in which a 10 2 cm" 3 peak density halo would collapse to 
form a star. As Fig. 8 shows, this would demand that such a 
dense minihalo be within tmbh ~ 240 pc of the halo hosting 



the first Pop III star at the time that this first Pop III star 
collapses to form a black hole. Such minihaloes would lie 
well within the H II region of the Pop III star and so may 
be strongly heated and partially ionized. However, the gas in 
such a halo is expected to remain neutral at densities n <; 2 
cm" 3 (Alvarez et al. 2006; see also Abel, Wise & Bryan 2006; 
Susa & Umemura 2006), which justifies the assumption that 
we use in our simulations of neutral neighboring minihaloes 
adjacent to the minihalo hosting the first Pop III star, for 
haloes containing gas at densities J; 10 cm" 3 . 



5 BLACK HOLE ACCRETION AND 
FEEDBACK 

In order for a black hole with a mass of ~ 100M Q , such as 
may have been formed from the direct collapse of a Pop III 
star at z ~ 20 (e.g. Heger et al. 2003), to attain a mass of ~ 
10 9 Mq, such as have been inferred for the black holes that 
power quasars at z i; 6, an average accretion rate of J> 1 M 
yr" 1 is required. Here we analytically estimate the accretion 
rate of a black hole that begins with a mass mi, the resulting 
mass of the black hole m as a function of time t since its 
formation, and the luminosity that would be produced in 
the process of this accretion as a function of time, for a 
given density and temperature of accreted gas. We assume 
that the black hole accretes gas from a cloud with uniform 
density and temperature and with dimensions much greater 
than the accretion radius of the black hole, r acc ~ Gm/c 2 , 
where c s is the speed of sound in the gas. In this case, we 
can estimate the Bondi-accretion rate as (Bondi 1952) 



m B 



2ir(Grn) 2 mH_n 



(19) 



Assuming that a fraction e ~ 0.1 of the accreted mass is 
converted into energy and radiated away, so that the black 
hole grows according to dm/dt ~ (1 — ejraa, integrating from 
t=0 and m = mi yields a time dependent mass of the black 
hole of 



2tt(1 - e)m^G 2 nt 
(3fcT)§ , 



(20) 



where T is the temperature of the gas at infinity and where 
we have used c s = 3fer/mH. 

Using equation (19) for the mass accretion rate, we com- 
pute the luminosity L generated in the process of the black 
hole accreting mass, operating as a miniquasar, as 



dm 



e dt 



2Tve(Gm) 2 mHnc 2 
(1 - e)c 3 



(21) 



where c is the speed of light. Expressing this as a fraction 
of the Eddington luminosity of the black hole yields 



L 



ecrTGmnc 
2(1 - e)c| 



(22) 



where <tt is the Thomson cross section. With equation (20) , 
we thus find 



to^Gm^nc 



1 2tt(1 - e)m|G 2 «i 



^Edd 2(l-e)(3Jfe B T)i \ m i 



(3fc B T)5 



.(23) 
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Figure 9. The mass of an initially 100 Mq black hole as a func- 
tion of time, assuming the black hole accretes gas at a tempera- 
ture of 200 K and at a constant density. The solid, dotted, dashed, 
and dotted-dashed lines show cases with densities of 0.1 cm -3 , 
1 cm -3 , 10 cm -3 , and 100 cm -3 , respectively. The triple dot- 
dashed line shows the mass of the black hole as a function of 
time, assuming that the black hole accretes at the Eddington 
limit. Clearly, for the black hole to begin accreting at near the 
Eddington limit while its mass is still of the order of 100 Mq, it 
must accrete gas that has a density > 10 2 cm -3 . 
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Figure 10. The luminosity emitted by a miniquasar fueled by 
accretion onto an initially 100 Mq black hole, expressed as a 
fraction of the Eddington luminosity of the black hole, as a func- 
tion of time. The densities assumed for the accreting gas are the 
same as those of Fig. 9, and the temperature is again assumed to 
be 200 K. 

To use our analytical argument to estimate the early 
growth rate of the first black holes, we need to determine 
the densities and temperatures of the accreting gas. We can 
derive their typical range of values from our simulations of 
the thermal and dynamical evolution of the gas in the vicin- 
ity of the newly formed MBH (in Sections 3 and 4) . We take 



it here that the accretion of DM particles contributes a neg- 
ligible amount to the overall mass accretion rate, making 
the assumption that DM particles are collisionless. Conse- 
quently, DM particles within the virialized halo in which the 
black hole resides cannot lose angular momentum efficiently 
enough to fall into the MBH. If, however, DM particles were 
self-interacting and were thus not collisionless, the accretion 
of dark matter could enhance the growth of the black hole 
(see Hu et al. 2006). 

We do not consider the effects of radiative feedback on 
the accretion rate of the black hole, although it has been 
shown that at high accretion rates the radiation emitted 
from quasars can inhibit further gas infall (see Di Matteo, 
Springel & Hernquist 2005; Springel, Di Matteo & Hernquist 
2005). Indeed, radiative feedback from the accretion of mat- 
ter onto the black hole may heat the accreting gas which 
surrounds the black hole to higher temperatures, which will 
in turn lower the accretion rate, as can be seen from equa- 
tion (19). This assumption is partially justified by the low 
accretion rates and associated luminosities encountered in 
our simulations, although it is dependent on the uncertain 
efficiency with which the accreting gas may absorb the emit- 
ted radiation. We have also assumed that the black hole is 
not moving with respect to the accreting gas, as might be 
the case if the black hole receives a kick at its birth. These 
last two assumptions make our results for the accretion rates 
onto the black hole upper limits. 

Finally, as for our assumption that the fraction of ac- 
creted mass that is converted to radiation, e = 0.1 will be 
an overestimate for accretion rates that are much below the 
Eddingtion accretion rate. In general, e may be much lower, 
and so our results for the luminosity generated by the accre- 
tion onto black holes should also be considered upper limits. 

Fig. 9 shows the mass of an initially 100 Mq black hole 
as a function of time, using equation (20), for various den- 
sities of the accreting gas and assuming a temperature for 
the accreting gas of 200 K, which is the temperature of the 
highest density gas we find in our simulations of a relic H II 
region merging with a pre-collapse un-ionized minihalo. Also 
shown in Fig. 9 is the mass of the black hole as a function 
of time, assuming that it accretes at the Eddington limit at 
all times. This figure clearly shows that in order for a black 
hole with a mass of the order of ~ 100 Mq to accrete at 
near the Eddington rate, it must accrete gas with a density 
<; 10 2 cm -3 . Since these densities are only found inside col- 
lapsing haloes, and are much higher than those of the gas 
in the expanding relic H II region around the black hole, for 
such a black hole to begin accreting at the Eddington rate 
it must merge with a minihalo that has a central density 
10 2 cm -3 . 

In the case of accretion at the Eddington limit, the black 
hole attains a mass of ~ 10 9 Mq in ~ 800 Myr, which means 
that such a black hole could marginally grow to become a 
SMBH by redshift z ~ 6 and power the observed quasars 
at that redshift if it forms from a collapsed Pop III star at 
a redshift of z <; 30 and accretes at or near the Eddington 
limit continuously. However, for densities lower than ~ 600 
cm~ 3 , or gas temperatures higher than ~ 200 K, an initially 
~ 100 Mq black hole likely cannot grow fast enough to ex- 
plain these observations, even if it attains mass also through 
mergers with other massive black holes (see Malbon 2006). 
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Black holes of initially higher masses could accrete at 
higher rates, since the accretion rate varies as m 2 , as shown 
in equation (19). An initially 500 M Q black hole would re- 
quire accreting gas densities J> 10 2 cm~ 3 at temperatures Ss 
200 K in order to accrete at the Eddington limit. Accreting 
at this limiting rate continuously, such a black hole could 
achieve a mass of 10 9 M within 700 Myr. Thus, such a 
black hole could be the progenitor of the SMBHs at redshift 
z ~ 6, if it merges with a dense minihalo and so begins ac- 
creting at or near the Eddington limit starting shortly after 
its formation by the collapse of a Pop III star at a redshift 
z £ 20. 

As Fig. 8 shows, mergers of sufficiently dense and suffi- 
ciently nearby pre-collapse minihaloes with a minihalo host- 
ing a Pop III remnant black hole may have led to accretion 
of gas with densities <; 10 2 cm -3 at temperatures of ~ 200 
K. Now, as Fig. 9 shows, if a 100 Mq remnant black hole 
accreted such high density gas during these mergers, it could 
have accreted the gas efficiently, the accretion rate even ap- 
proaching the Eddington limit. However, these high gas den- 
sities are confined to only the regions very near the centers 
of the merging minihaloes, as can be seen in Fig. 5. Indeed, 
for the case of the merger shown in this figure, in which the 
peak density at the center of the merging minihalo becomes 
n ~ 10 3 cm~ 3 after 100 Myr, the remnant black hole ac- 
cretes only much lower density gas, residing in a region with 
a gas density of only ~ 10~ 2,5 cm~ 3 at this same time in 
the course of the merger. The accretion history of a Pop III 
remnant black hole thus sensitively depends on the location 
of that black hole, as well as on the highest densities of the 
gas involved in the merger. 

One generic requirement for a Pop III remnant black 
hole to accrete at the Eddington limit early enough to 
achieve a mass of ~ 10 9 Mq by z ~ 6 appears to be that it 
accrete the high density gas within a merging pre-collapse 
minihalo which can lie only a few hundred parsecs from the 
black hole at the time of its formation, as Fig. 8 shows. 
Another requirement is that the peak density in this pre- 
collapse halo must be at least n ~ 10 cm~ 3 at the time 
when the original Pop III star stops emitting LW photons 
and collapses to become a black hole, so that the density 
in this halo can rise to 10 2 cm~ 3 . As Fig. 9 shows, this 
is the density that the accreted gas must have in order for 
accretion to take place at the Eddington limit for a black 
hole with a mass of the order of 100 Mq. If either of these 
requirements is not satisfied, and the black hole does not 
merge with a sufficiently dense or nearby minihalo, then the 
black hole will not begin accreting at the Eddington rate 
before it is incorporated into a larger halo that can gravi- 
tationally bind the hot gas in the relic H II region, and so 
begin to raise the density of this gas again. The time scale 
for this to occur is ~ 100 Myr (see Yoshida 2006), which 
means that there would be a substantial delay of the order 
of the Hubble time at a redshift of z ~ 30 before a Pop III 
relic black hole could begin efficiently accreting, unless it 
merges with a dense minihalo soon after its formation. 

Using equation (23), we find time-dependent luminosi- 
ties of a miniquasar fueled by a black hole remnant of a 
Pop III star, expressed as a fraction of its Eddington lu- 
minosity, as shown in Fig. 10. Here we have assumed the 
same range of values for the density, and the same tempera- 
ture, of the accreted gas as was assumed in Fig. 9. The same 



requirements for the black hole to begin accreting at the Ed- 
dington limit apply to the miniquasar, fueled by accretion 
onto the black hole, to radiate at the Eddington luminos- 
ity. Thus, the Eddington luminosity will only be achieved 
within ~ 100 Myr for the case of a merger with a neigh- 
boring pre-collapse halo with a central density above ~ 10 
cm~ 3 , and which lies within a few hundred parsecs from the 
black hole, at the time of the formation of the black hole. 
As Fig. 10 shows, black holes that do not accrete cold gas 
at densities <; 10 cm will not radiate efficiently for well 
over a Hubble time at redshifts z <; 20. Therefore, it may be 
difficult for miniquasars to be strong ionizing sources which 
contribute substantially to reionization at redshifts J> 15, 
unless they are fueled by the accretion of dense gas onto 
Pop III seed black holes that merge with dense minihaloes 
soon after their formation (see Madau et al. 2004; Ricotti 
& Ostriker 2004; Kuhlen & Madau 2005). We note, how- 
ever, that the number density of accreting black holes at 
these high redshifts is unknown, and if this number density 
is high enough then even inefficiently radiating miniquasars 
may be able to substantially reionize the Universe. 



6 SUMMARY AND CONCLUSIONS 

We have carried out three-dimensional numerical simula- 
tions of the evolution of the first relic H II regions and of the 
mergers of pre-collapse minihaloes with the minihaloes that 
host the first Pop III remnant black holes at the centers of 
the first relic H II regions. Although the first Pop III stars 
which form in minihaloes emit enough ionizing radiation to 
evacuate the primordial gas that resides in these haloes and 
to disperse it into the IGM, we find that the radiative feed- 
back from the first stars, nonetheless, does have important 
effects that could help to foster continued star formation in- 
side the relic H II regions that they leave behind. Namely, 
the recombination of these relic H II regions, due to their 
dynamical expansion, does not proceed to completion for 
over a Hubble time at the redshifts at which the first stars 
likely formed, z i> 20. The residual electron fraction that 
persists for long after the death of the central star in these 
regions allows for continued molecule formation, and thus 
high fractions of both H2 and HD are produced. 

Due to the large size of the first H II regions, of the 
order of a few kpc (e.g. Alvarez et al. 2006), and to the high 
H2 fraction that is generated in the recombination of the 
ionized primordial gas after the death of the central star, 
of the order of 10~ 3 , we find that shielding of molecules in- 
side the first relic H II regions could have been significant, 
providing an optical depth to LW photons of the order of 
10 (see also Haiman et al. 2000; Ricotti, Gnedin & Shull 
2001; Oh & Haiman 2002). With the effects of continued H 2 
formation and of high optical depth to LW photons taken 
together, we have shown that molecules inside the first relic 
H II regions could have persisted at least down to redshifts 
of z ~ 15 before the background LW flux in the Universe 
grew strong enough to efficiently photodissociate H2 in these 
regions. We have further shown that minihaloes with cen- 
tral densities of ii 10 2 cm~ 3 that suffer complete destruc- 
tion of H2 molecules due to photodissociating flux from a 
Pop III star can reform H2 and collapse within 60 Myr af- 
ter the death of the Pop III star (see also Susa & Umemura 
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2006; Mesinger et al. 2006). Thus, second-generation Pop III 
stars could have formed inside relic H II regions within mini- 
haloes with central densities of at least this order of magni- 
tude, and, furthermore, these sites of second-generation star 
formation could have been shielded from the background 
molecule-dissociating LW flux established by star formation 
elsewhere in the Universe. 

In addition, we have shown that a very high fraction 
of HD molecules forms inside the first relic H II regions, 
reaching levels of the order of 10~ 7 . As this fraction is well 
above the critical fraction of HD needed for efficient radia- 
tive cooling of the primordial gas to the temperature of the 
CMB, when this gas is in LTE, in principle, the first relic 
H II regions produce sufficient HD to allow the formation 
of Pop II. 5 stars. However, we emphasize that the primor- 
dial gas in the first relic H II regions must become gravita- 
tionally bound inside larger mass haloes in order for it to 
become dense enough to collapse and form stars, if indeed 
Pop II. 5 stars are to form inside these regions. Nagakura 
& Omukai (2005) investigate this possibility and find that 
stars could have formed from gas cooled by HD to tempera- 
tures ^ 100 K inside haloes massive enough to bind dynam- 
ically expanding relic H II region gas. Three-dimensional 
numerical simulations of such larger mass systems are nec- 
essary to more reliably discern whether Pop II. 5 star for- 
mation can occur in the relic H II regions surrounding the 
first stars. Interestingly, recent three-dimensional radiation- 
hydrodynamics simulations reported by Yoshida (2006) sug- 
gest that the fragmentation scale of HD-enriched gas inside 
relic H II regions can become of the order of 10 Mq, which 
corresponds to the characteristic mass expected for Pop II. 5 
stars (Johnson & Bromm 2006). 

The origin of the first SMBHs, which are observed to 
have masses of the order of 10 9 Mq at z ~ 6, is an impor- 
tant and challenging question to address. In order for Pop III 
remnant stellar black holes, with initial masses of the order 
of 100 Mq to grow to become the observed SMBHs, they 
must begin accreting at or near the Eddington rate very 
soon after their formation. Indeed, in order to explain the 
growth of supermassive black holes, much recent work draws 
on the assumption that this is the case (e.g. Malbon et al. 
2006; Li et al. 2006). In particluar, 100 and 500 Mq black 
holes must accrete at the Eddington limit for 800 and 700 
Myr, respectively, in order to grow to 10 9 Mq. Addition- 
ally, recent work by Malbon et al. (2006) shows that it is 
indeed the accretion of gas that largely fuels the growth of 
massive black holes at redshifts z *J 2, although at lower 
redshifts black hole mergers can become the dominant pro- 
cess for black hole growth. Because the radiation from the 
progenitor Pop III stars which collapse to form such black 
holes evacuates the gas from their host haloes, mergers with 
neighboring minihaloes must occur in order for the remnant 
black holes to accrete high enough density gas to accrete at 
or near the Eddington limit. 

Our results indicate that there may have been a sub- 
stantial time delay between MBH formation and the onset 
of efficient accretion in situations where Pop III seed black 
holes do not merge with neighboring neutral minihaloes with 
densities that are J; 10 cm~ 3 at the time of black hole for- 
mation. In this case, efficient accretion cannot begin until 
the black hole is incorporated into a halo sufficiently mas- 
sive to gravitationally bind the hot relic H II region gas (see 



Yoshida 2006), and this can take of order the Hubble time at 
z ~ 30, or ~ 100 Myr. It remains therefore an open question 
whether the SMBHs that power the z J> 6 quasars could have 
grown from stellar remnants, or whether more massive seed 
black holes were required (e.g. Bromm & Loeb 2003; Lodato 
& Natarajan 2006). More realistic, fully cosmological simu- 
lations which accurately track the dynamics and growth of 
accreting black holes over ~ 10 9 yr will be required to fur- 
ther clarify the possible connection to the SMBHs observed 
at lower redshifts. 

If the conditions we have found for accretion at the 
Eddington limit are not met, then the miniquasars fueled 
by accretion onto the Pop III remnant black holes will also 
not emit radiation at the Eddington limit until well after 
the formation of the black holes. This possibility of ineffi- 
cient accretion onto these black holes has important conse- 
quences for theories of the reionization of the Universe by 
miniquasars, as it is often assumed that accretion at the Ed- 
dington rate is achievable for miniquasars at redshifts z *J 
15 (see Madau et al. 2004; Ricotti & Ostriker 2004; Kuhlen 
& Madau 2005). Although, a sufficiently high number den- 
sity of accreting black holes could still lead to significant 
reionization by miniquasars. If accretion onto miniquasars 
is inefficient, as discussed in Madau & Silk (2005), the gen- 
eration of the observed near-infrared background excess by 
Pop III stars (e.g. Santos, Bromm & Kamionkowski 2002; 
Kashlinsky 2005) may not necessarily imply an overproduc- 
tion of X-rays emitted by miniquasars fueld by accretion 
onto Pop III remnant black holes. However, copious X-ray 
production could result from accretion of gas from a bi- 
nary companion (see Belczynski, Bulik & Rudak 2004; Saigo, 
Matsumoto & Umemura 2004), even if the binary system re- 
sides in a low density environment. Future work will have to 
include more detailed, fully cosmological simulations with 
a prescription for radiative feedback from Pop III remnant 
black holes to self-consistently address these issues. 
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